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Abstract 

We derive the pressure tensor and the heat flux to accompany the new macroscopic 
conservation equations that we developed in a volume-based kinetic framework for 
gas flows [1]. This kinetic description allows for expansion or compression of a fluid 
element, and is characterized by a flux of volume that conventional modelling does 
not account for. A new convective form of transport arises from this new approach, 
which is due purely to macroscopic expansion or compression. It supplements the 
conventional transport processes, which are transport due purely to mass motion 
(classical convective transport) and transport due purely to random motion (dif- 
fusive transport). In addition, we show that the thermodynamic parameters of the 
fluid (temperature and pressure) appear with new non-equilibrium contributions 
that depend on density variations in the gas. Our new model may be useful for de- 
scribing gas flows that display non-local-thermodynamic-equilibrium (rarefled gas 
flows), flows with relatively large variations of macroscopic properties, and/or highly 
compressible fluids and flows. 

Key words: gas kinetic theory, Boltzmann equation, compressible fluids and flows, 
non-isothermal flows, microflows, Navier-Stokes equations, rarefled gas dynamics, 
constitutive relations 



1 Introduction 



The conventional fluid mechanics description is the set of three Navier-Stokes 
equations, which can be retrieved from the kinetic theory of rarefied gases 
using, for example, a perturbation method (known as the Chapman- Enskog 
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expansion) truncated at first order p[3] . The validity of the Navier-Stokes set 
of equations is predicated on small deviations from thermodynamic equilib- 
rium, which can be interpreted physically as relatively slow and low amplitude 
variations of macroscopic variables within the fiows. Describing fiows beyond 
the broad range of applicability of the Navier-Stokes model (such as hypersonic 
fiows, or micro or nano scale gas fiows) remains an active area of investigation 



An intuitive way of constructing hydrodynamic models beyond the range of ap- 
plicability of Navier-Stokes is high-order development of the Chapman-Enskog 
expansion. These are the Burnet and Super-Burnet equations. However, these 
high-order expansions have various and well-known physical inconsistencies, 
including frame- dependent diffusive fluxes [9j, violation of the second law of 
thermodynamics [TU], and other instabilities [TT]. Various ways of overcoming 
these inconsistencies have been proposed [T^|IT3] . 

In a previous paper [Ij, we discussed the representation of fluids and flows 
from both the gas kinetic and the classical fluid mechanical viewpoints. We 
introduced a kinetic representation of dilute gases based on volume variations. 
The fundamental change introduced in our new model was that, while gases are 
physically represented by particles undergoing random motions and collisions, 
the free volume around each gaseous molecule should be taken into account 
in their microscopic modeling, and that this free microscopic volume around 
each particle is linked to macroscopic variables such as the density of the 
fluid. We modifled Boltzmann's original kinetic equation, which is a simple 
balance of the number of particles, to account for the microscopic free volume 
around particles. We therefore arrived at a new set of macroscopic conservation 
equations. In addition, two different macroscopic velocities, leading to two 
different definitions of kinetic peculiar velocities, were distinguished: the mass- 
velocity (which is the velocity of the centre of mass of an element of fiuid), and 
the volume- velocity (which allows for expansion/compression in the motion of 
a fiuid element). 

The aim of the present paper is to build on this new volume-based description 
by developing first-order constitutive relations for the pressure tensor and heat 
fiux to accompany it. 

This paper is structured as follows: 

Section [2] summarises our new volume-based kinetic approach. 

In section [31 an unconventional convective form of transport arising from our 
new kinetic approach is presented. This convective transport is due purely to 
macroscopic expansion/compression, and supplements the two conventional 
transport processes, namely a transport due purely to mass motion (the clas- 
sical convective transport) and a transport due to random motion (without 
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any bulk motions, and which is called diffusive transport). 

In section H] we derive new constitutive relations. This derivation is mainly 
based on the identification of diffusive fluxes as convective fluxes in our new 
framework. We argue that actual diffusive fluxes are those due purely to parti- 
cle random motions, so they should be fluxes from which not only mass motion 
but also macroscopic expansion/compression motion contributions should be 
removed. A complete new set of hydrodynamic equations is then proposed. 

In section m we analyse the new quantities involved in our proposed hydrody- 
namic model. We show that the flux of volume, which is normally a convec- 
tive flux, behaves also like a diffusive flux. Finally, the new thermodynamic 
parameters generated by our model are compared with the thermodynamic 
parameters of the classical description. 

(Notation: In this paper, the square of the modulus of a vector Z we denote 
by 



2 A volume-based kinetic description of gas flows 

In the kinetic theory of simple monatomic gases, the gaseous particles are 
only characterized (at any time t) by their position, X, and velocity, ^. In 
order to describe situations sensitive to variations of macroscopic parame- 
ters, especially density, we introduced the microscopic free volume, f , around 
each gaseous particle into the set of microscopic parameters characterizing the 
molecules [I]. A distribution function f(t,X,^,v) is then a probability num- 
ber density of particles which, at a given time t, are located in the vicinity 
of position X, have their velocities in the vicinity of velocity ^, and have an 
assigned microscopic volume in the vicinity of volume v. 

As the dynamics of the gaseous particles is described through their point- 
masses, the kinetic equation describing the evolution of f{t,X,^,v) is the 
following Boltzmann-like equation: 

§|- + (e ■ V)/ + (F ■ V^)/ + =J J irf* - ffiMrd^d^,, (1) 

in which the term on the right-hand-side is the classical Boltzmann collision 
integral; / = f(t,X,C,,v) and /i = f(t,X,C,i,vi) refer to post-colhsion parti- 
cles, /* = f{t, X, V*) and fl = f(t, X, v^) refer to pre-coUision particles, 
= ^ — ^1 is the particle relative velocity, a the collision differential cross 
section, an element of solid angle, and denotes the formal operator 
= d/d^x + d/d^y + d/d^z- On the left-hand-side appears a new term in- 
volving W , which arises from the introduction of the new variable v into the 
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distribution function. This term represents the internal rate of change of vol- 
ume occupied by the gas [1]. In the following we will suppose no body force, 
i.e. F = (although situations of non-vanishing body force can be easily 
integrated into the description). 



2. 1 The macroscopic parameters of the gas 



The macroscopic parameters according to our modified distribution function 
f(t,X,$,,v) are defined as follows: 

• The local number- density of the molecules is given by 

/ f{t,X,^,v)d,d^. (2) 

-oo Jo 

This number density refers to the fixed observer's reference frame with co- 
ordinate elements (Xi,X2,X3), in which any motion and dynamics are de- 
scribed. 

• The local mean value, Q{t,X), of any property Q can be defined by, 

/+00 r-\-co 
/ Qf{t,X,^,v)d,d^ . (3) 
-oo Jo 

For example, the local mean microscopic volume, v{t, X), around each gaseous 
particle is given by. 



n(t, X)v{t, X)= / vfit, X, v)d,d^ . (4) 

J-oo Jo 

• From this mean value of the microscopic volume, a local mean value of the 
mass-density in the vicinity of position X is defined through: 

, n(t,X)M M 

n{t,X)v[t,X) v[t,X) 

where M is the molecular mass. 

It may be noted that in gas flows where the molecules can be assumed always 
uniformly dispersed (locally), the mean value of the microscopic volumes, v, 
is simply given by dx/{ndx) = 1/"^, regardless of definition (jl]). In this case 
the definition of density, equation ([5]), becomes p = Mn. However, we are 
interested here in general situations where these restrictions on the volume 
occupied by a group of particles are relaxed; this is in order to get more 
efficient descriptions of flows with density variations. While the density Mn 
can be regarded as the density of a gas without density variations, the actual 
mass-density of the general gas (i.e. the mass divided by the actual volume 
occupied by the gas) is given by p in equation 
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Two mean velocities can be defined. First, the local mean mass- velocity, 
Um(t,X), is given through 

Mn{t,X)UUt,X) = J J MUXt,X,^,v)d^d,. (6) 

A local mean volume- velocity, U^(t,X), can also be defined: 

v{t, X)n{t, X)U,{t, X)=J J v^fit, X, v)d^d,. (7) 

While Um{t, X) is purely the velocity of the centre of mass of a "fiuid particle" 
(or an element of fiuid) around X, the velocity Uy{t,X) takes into account 
any expansion or compression velocity of this fiuid particle. This is because 
the weighting element in equation ([7]), which is a volume, is a microscopically 
changing property in contrast to the mass in equation ([6]). 

From these two local mean velocities, two peculiar velocities can therefore be 
introduced. The usual mass-velocity definition of peculiar velocity is 

C = ^-Um. (8) 

But another peculiar velocity may be given through the volume-velocity, i.e. 

C' = ^-U,. (9) 



2.2 A new set of macroscopic conservation equations 



Multiplying the modified kinetic equation ([T]) by the four microscopic quanti- 
ties V, M, (M^), (M^^) and then integrating over d^ and d^, we obtained the 
following set of macroscopic conservation equations: 



Continuity 
Mass-density 

Momentum 
Energy 



Dn 

lyt 



-nV ■ Ur, 



Dt MVn ^ 



W 



DUr, 



n- 



Dt 



n- 



D 
Dt 



-U + e- 



VP, 



-V ■ iP ■ u„ 



V ■ q 



(10) 

(11) 

(12) 
(13) 



where the material derivative D/ Dt = d/dt + Um ■ V. 
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In this set of equations J„, P, Cm and q simply denote the following quantities: 

Jv = [ f vCfd^d^ , (14) 



P 



ne. 



J J QC,fd,d^ , (15) 

rn = lj J C'fd^d, , (16) 

q=^/ Jc'Cfd^d,. (17) 

It is important to note that these quantities, which are all written with the 
peculiar velocity C based on the mass- velocity Um, are simply quantities ap- 
pearing through the rational mathematical derivation of the macroscopic equa- 
tions. They are not yet shown to be the same as the pressure tensor, heat flux, 
or internal energy in the fluid. In particular, the presence here of C in equations 
(1151) to f|T71) is due to the definition within equation ([6]) and to the Boltzmann 
integral collision which vanishes when (M^) and (M^^) are introduced in the 
integral forms of equation ([T]). 

Our new set of macroscopic conservation equations contains four equations 
instead of the usual three. Equation (fTO!) simply ensures the conservation of 
mass. It involves a number- density, n, which should be regarded as accounting 
only for the number of particles, and the mass- velocity, Um- On the other hand, 
equation (ITT]) gives the evolution of the actual density p as a thermodynamic 
parameter of the fluid which varies with the volume occupied by the fluid. 



3 An unusual convective transport driven by a flux of volume 



Using the definitions of the mean velocities and peculiar velocities, equations 
(E]) to ([9]), the following identities can be easily proved: 

C'fd^d, = -n{U, -UJ , J J vC'fd^d, = , (18) 

and 

nvUy = J J vUmfd^d^ ^ J J ^^f'^i'^-" ■ 
So equation f|T^ may be re-written: 

3^ = nv{U,-Um) . (20) 

In equation (l20l) . J^, is characterizing a macroscopic motion defined by ([/„ — 
Um), given that Um and [/„ are both macroscopic velocities by definition. The 
flux therefore acts like a convective flux in which the convected element is 
the volume transported at the macroscopic velocity {U^ — Um)- 
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Fig. 1. Schematic of fluid element mass velocity, Um, and volume velocity, 17^; a) 
the macroscopic motion of a fluid element in a fluid without density variations, b) 
macroscopic motion in a fluid with density variations, C/„ = Um + l/{nv)Jy. Each 
of the volumes contain the same number of molecules. 



This unconventional macroscopic motion represented by J^, corresponds to a 
macroscopic change of volume of a "fluid particle" in hydrodynamic descrip- 
tions. A fluid particle is not a simple point-mass, but a volume element; this 
is illustrated in Figure [H If 6V represents the macroscopic change of volume 
of a fluid element over time 6t, then an approximate expansion/compression 
speed is given by, || {nv)''^J^ ||^ 5{V^^^)/5t , where || . || denotes the modulus 
of a vector. The velocity J„ is oriented in the direction from high density to 
low density in the case of an expansion, and vice versa for compression. 

The macroscopic motion characterized by the flux J„ affects the classical in- 
terpretation of "convective fluxes" . In our volume-based description, the term 
"convective flux" should include the macroscopic motion deflned through the 
flux J^, as well as the traditional centre of mass velocity Um- 



4 Constitutive relations in our volume-based description 

In this section we identify the various fluxes to be associated with the pressure 
tensor and heat flux, and the expression of internal heat energy. Then we derive 
constitutive relations for these various quantities, and for the flux J^. 
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4-1 An expression for the flux of volume 



Let us regard the microscopic volume, v, as a property attached to the molecules 
in a gas that has no macroscopic motion (or macroscopic expansion or com- 
pression). Then equation fll4p . as it stands, represents a diffusive ffux in which 
the microscopic element undergoing diffusion is v. By a diffusive ffux we un- 
derstand a ffux due purely to the random motions of the particle, without any 
form of bulk; motion. In this case, following phenomenological laws of diffusion 
we may expect the diffusive volume ffux to have the form 



= —Ky Vf , (21) 

or, equivalently, by using equation ([5]), 

j^ = ^S/p, (22) 
P 

where we will call k„ the "volume diffusion coefficient" . 

On the other hand, we showed in the previous section that the ffux phys- 
ically characterizes the macroscopic motion displayed in equation (l20l) . Any 
ffux which is given by this macroscopic motion should be taken as a convective 
ffux. So, the diffusive-like expression for J^, in equation (pT]) may be interpreted 
as a "fictitious" diffusion of volume. However, we resolve this dual nature of 
Jy by explicitly proving in section Oof this paper that, in fact, in addition to 
defining a new convective transport, Jy corresponds to a particular diffusion 
of mass. That it has both a convective and a diffusive nature justifies equation 
(I22D. 



4-2 The pressure tensor, heat flux and internal energy 



Classical kinetic theory attributes the pressure tensor and heat ffux to ffuxes 
deriving from the peculiar velocity p|3j . It is also generally understood that 
any macroscopic motion is defined only by the mass- motion, Um- So purely 
random components of the partitioning motion would then be justifiably de- 
fined by the peculiar velocity C through equation ([8]). The concept of pure 
macroscopic expansion/compression motions, as pointed out in section [31 does 
not exist. 

In the conceptual framework of hydrodynamics (and kinetic theory), a diffu- 
sive ffux should be a purely internal ffux of momentum or energy in which 
contributions due to any macroscopic motions are removed. This requirement 
is not satisfied by the three quantities defined in equations (fT5|) to (fTTll . The 
macroscopic motion that arises from equation (l2Up . and then the convective 
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contributions due to J^, are still inherent in fluxes expressed using the peculiar 
velocity C . 

However, all macroscopic motions — including macroscopic expansion and 
compression of a fluid element — are removed from the peculiar velocity C 
deflned in equation ^ because this equation can be re-written, using equation 
([20]), as 

C = i-Um- — J. . (23) 
nv 

Furthermore, equations ([8]), ([9]) and fl20|) yield 

C -C = U,-U^ = — J, . (24) 

nv 

We could therefore express the fluxes q and P, and the quantity from 
equations (ITSl) to (fTTI) . which are simply quantities deflned with the peculiar 
velocity C, using the following decomposition: 

C = C' + {C -C') = C' + {U, - Urn) . (25) 



Then P{t,X) from equation (JTSi) with equation (l25l) becomes 

Pij — P ij TTTT ( Ji) J-u)i7 5 (^6) 

where J^J„ = (J^J„)ij is the second order tensor obtained by the product of 
the coordinate components of vector J^, i.e. {3^Jy)ij = (J^)j(J^,)j, and P' is 
the actual momentum flux tensor 

P% (t,X)= J jc[C'Jd,d^. (27) 

A scalar pressure p' can then be introduced through the sum of the three 
diagonal terms of P': 

3p' = MV'ii . (28) 



From equations (fT6l) and (!25l) ej„ becomes 

1 



2(nt;)2 " ' 

where the actual internal energy is given as 



(29) 



^-J J C'fd^d, . (30) 



A temperature T' can also then be deflned: 



hr = Me[^ , (31) 
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and using equations (128|) . (l30l) and (13T!) we obtain 

p' = nfcT' . (32) 

Using equation fl2^ . the flux q, defined in equation f[T7l) . is re- written 

q{t,X) = ^J J[C' + {U,-Um)fC'fd^d, + {U,-Um)^J Jc'fd^d,. (33) 
Then equations (HHD, (El) and (I29l) yield 

q = q' + — P' ■ J„ - ^T^JP. + - ( - 7^7^ J'j J. , (34) 
nv 2v [nv)^ v \ 2[nvj^ J 

which can be written 

q = q' + — P'- J^ + - (e-„--^J^) J„ , (35) 
nv v \ (nf J 

where the actual energy flux vector, q' = q'(t,X), is given by 

q' = ^ / / C'C'fd^d, . (36) 

The fluxes q' and P', and the quantity e[^, appear as new internal quantities 
with contributions due to any bulk motions removed. The "actual" diffusive 
fluxes, by which we mean those with any kind of convective fluxes removed 
(including those described in section [3]), are the new quantities P' for the 
momentum and q' for the energy. In other words, diffusive fluxes in our volume- 
based kinetic description are those expressed using the peculiar velocity C, 
in accordance with equation fl23|) . and not the fluxes appearing in equations 

dE]) to ([HD . 

Equation (!29l) describes the total internal energy of a fluid particle, Cm, as 
being composed of a "heat energy", e-„, and an elastic energy component 
driven by J^,. If we imagine a fluid particle (containing a certain amount of 
mass) moving from a high density region to a low density region in the fluid, 
this will cause a dilatation of the fluid particle if the quantity of mass in the 
fluid particle is not allowed to change. Particles may move out of or into a 
fluid particle, but the amount of mass should remain the same, cf. Figure [TJ 
According to the first law of thermodynamics, the fluid particle heat energy 
must exclude any kind of potential energy. In equation (!29|) the term in is 
a kind of potential energy because it is related to the density gradient, and so 
represents the dilatation work done by the fluid particle in the fluid by density 
variations. So the actual heat energy should correspond to e-„, not Cj^. 
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The process we have just described does not exist within classical fluid me- 
chanics, where the total internal energy of the fluid particle is associated only 
with the heat energy which is directly related to the temperature. For a moving 
and expanding or contracting fluid particle, the classical fluid mechanics defi- 
nition of internal energy does not agree with the first law of thermodynamics 
conception of heat energy as it combines a potential energy of the fluid particle 
with its actual heat energy. We show below that our distinction of the heat 
energy of the fluid may become sensitive in some extreme flow configurations 
(e.g. rarefied gas flows under high temperature or density gradients). 



4.3 Constitutive model and the final set of hydrodynamic equations 



As the fluxes due purely to random molecular motions in our kinetic volume- 
based description correspond to q' and P', it follows from phenomenological 
laws of diffusion that these fluxes may be modelled at first order by linear 
constitutive relations. These relations are the classical Newtonian law for the 
stress tensor within the pressure tensor, P', and Fourier's law for the heat 
flux, q'. The corresponding temperature and pressure in the fluid are then T' 
and p', respectively. 

As a result, we obtain a set of hydrodynamic equations given by the four 
macroscopic equations closed by equations fl26l) . (!29|) and (!35|l . in which: 



MP,=p*,,-na^+a^j-',g^%. (37) 




with yu' a dynamic viscosity, k'^ a heat conductivity, rj' a bulk viscosity, a 
volume diffusion coefficient, all to be determined in this volume-based frame- 
work. 

The macroscopic motion which defines the momentum flux P' as a diffusive 
flux is Uy. It follows that the phenomenological law of diffusion used to express 
P' should be applied with a gradient taken over U^, as we write in equations 
( 1371) . So momentum diffusion is not generated by the gradient of the mass- 
velocity Urn alone, but by the dilatation motion (f/„ — Um)- 

In the Appendix to this paper, we comment further on the expression of the 
stress-tensor in equations (1371) and compare it to the common fluid mechanics 
stress tensor that involves only the mass- velocity Um- 
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The complete set of new hydro dynamic equations is rewritten below, for con- 
venience: 



Continuity 



Dn 
15t 



-nV ■ U„ 



Mass-density 



Dp 
Dt 



n 



M 

Momentum 

DU^ 

n = 

Dt 

Energy 
D 



n 



V ■ [J.] - W 



-V- P' 



nv 



2 J. J. 



n- 



Dt 



2^m 




(3J 



(39) 



(40) 



(41) 



These equations are solved for five unknown fluid macroscopic parameters: 
the number-density, n, the mass-density (which involves the actual volume 
occupied by the gaseous particles), p, the mass- velocity, Um, the pressure, 
p', and the internal energy, (or, equivalently, the temperature, T'). The 
equation set is then closed by the constitutive model equations (137|) . and 
Me^„ = (3/2) fcT' , p' = nkV and U^ = U^ + (nt;)-iJ„. 



It is important to note that in the momentum equation PUI) . the macroscopic 
velocity on the left-hand-side is the mass- velocity, Um-, while on the right-hand- 
side, according to equations (|57l) . appears the volume velocity, U^. The modi- 
fications introduced by our new hydrodynamic description, when compared to 
the usual Navier-Stokes set of equations, recall the recent phenomenological 
modifications of Brenner [15]. That is to say, replacement of the mass velocity 
by the volume velocity in (only) the stress tensor, and the contribution to the 
heat flux of the work done by any variation of volume. 

In addition to Brenner's original modifications, second order terms due to the 
volume flux appear in both the momentum and the energy equations. Another 
difference between our model and Brenner's is our distinction between the 
real mass-density of the fluid, which we denote p, and the mass-density of an 
equivalent incompressible fluid, which is Mn in our notation. 
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In contrast to classical kinetic theory, equation (132|) is no longer an equation 
of state, as n does not encompass the real volume of an element of fluid. This 
equation will become a perfect gas equation of state only in the equilibrium 
situation, where Mn becomes p, i.e. if no density variation occurs. Variations 
of the density with the fluid macroscopic parameters are described through 
equation fl39l) . which is therefore functioning like an equation of state. 

An expression suggested for W in our previous paper [Tj, leads to the following 
re- writing of the density equation fl39|) : 

i§ = -f V.|J.l-« + ■ VT') +x + ■ Vp') , (42) 

p Dt Mn \ Ut nv J \ Ut nv J 

in which we may take a ^ 1/T' and x ~ K can be seen that when the flux 
Jd vanishes, e.g. in the case of a (local) thermodynamic equilibrium, equation 
fH2l) is of the form of the classical perfect gas equation of state. However, 
strictly speaking, the type of equation of state which should be displayed by 
the density equation for a vanishing flux of volume depends on how the rate 
of change of volume, W, is described in the kinetic description [T]. 

In equations fl38l) . fHOl) and fHTj) . the presence of the number density n (or the 
apparently constant density Mn) on the left-hand-side, is founded by a "mass- 
based integration" during the derivation of the conservation equations [1]. This 
is interpreted as following a fixed amount of mass when writing conservations 
of mass, momentum and energy. This is consistent with the material derivative 
appearing in all the left-hand-sides of these equations. The evolution of the 
domain occupied by the fixed amount of mass is contained in an equation of 
volume — that is our equation of mass-density that functions like an equation 
of state. 

Finally, it can be easily seen that our new hydrodynamic model becomes the 
Navier-Stokes set of equations when there is no mass-density variation, i.e. 
when J,, = 0. 



5 Analysis of the new hydrodynamic parameters 

5.1 The flux of volume as a diffusion of the mass-density 

Here we show that refers also to a particular mass diffusion process that 
justifies the expression for it which we hypothesised in equation (!22|) . 

In our new volume-based kinetic description, all diffusive fluxes have any 
bulk motion removed, including the expansion/compression motion defined 
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by (Uv — Urn). These diffusive ffuxes are therefore associated with the pecuhar 
velocity C given in equation (!23|l . Let us consider the following diffusive flux, 



J J pC'fd^d^ = lJ J MC'fd^d^ . (43) 



According to the definition of diffusive fluxes, the flux Jp corresponds to a 
diffusive flux of mass; strictly speaking, the element undergoing diffusion in 
this equation is the mass-density p. 

Using the definitions of the mean velocities, equations ([6]) and ([7j), with the 
definition of the peculiar velocity C in equation it can be shown that 

Jp = -pn{U, - . (44) 

On the other hand, from equation fl2^ we have 



U,-Um = —J. . (45) 
nv 

Combining equations (jH]) and fH51) . we obtain 

-Jp- = - -J. . (46) 
p V 

Equation (H6|) is a symmetrical relation between the flux of volume, J^, and 
the diffusive flux J p. 

As Jp is a diffusive flux, according to equation (H3l) . we can apply the phe- 
nomenological laws of diffusion: 

Jp = -fi^mVp , (47) 

where is the mass diffusion coefficient (strictly speaking, the mass-density 
diffusion coefficient). From equation fH6l) we then obtain, 

J, = ^Vp , (48) 
P 

which is the same as the expression for J„ written in equation fl22|) . but with 
the volume diffusion coefficient equal to the mass diffusion coefficient. 



In equation fl43|) or fl47j) . we see that the property of the fluid concerned in 
this diffusion is strictly the mass-density. In other words, this diffusion process 
acts to smooth out density inhomogeneities in a flow. This diffusion process 
should not be confused with self-diffusion, which is normally considered as 
a diffusion of matter through itself (in equilibrium). Self-diffusion occurs in 
homogeneous matter, with zero mass-density gradients, where one can imagine 
certain marked particles moving through similar but unmarked particles [3]. 
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Considering the mass-density as a "concentration of mass within the volume" , 
i.e. a mass divided by the volume of fluid, this diffusion process could be 
describes as "diffusion of mass through the volume of fluid". Or, from the 
symmetrical relation fl46p . "diffusion of the fluid volume through the mass". If 
we consider a gas at rest in which a density gradient is maintained, then this 
diffusion process may be regarded as a constant flux of volume imposed in the 
domain, or an apparent movement of the mass. On the other hand, if imposed 
density gradient is suddenly relaxed, then this diffusion corresponds to actual 
movement of mass from the high mass-density regions to low mass-density 
regions. 

Finally, regarding an element of fluid, any diffusive flux is measured in a 
reference frame that moves not only with the mass velocity but also undergoes 
expansion or compression with the element of fluid. In this global reference 
frame, the expansion or compression motion itself always looks like a diffusion 
of the mass-density. This is consistent with the expansion or compression 
associated with the flux J^, in earlier sections of this paper. The coefficients 
Km and Ky, which have the dimension of (length)^/time, resemble expansion 
or compression coefficients, rather than purely diffusion coefficients. 

5.2 Non- equilibrium thermodynamic parameters 

We now show that the scalars T' and p' defined in our volume-based description 
correspond to the gas thermodynamic parameters in non-equilibrium condi- 
tions. They coincide with the conventional scalars, T and p, only in equilibrium 
conditions. 

We introduce temporarily a scalar, T, which is related to the total internal 
energy, etn, in a similar way to equation (1311) . and a scalar p related to the 
sum of the three diagonal terms of the total momentum flux tensor P, in a 
similar way to equation (I28p . i.e., 

hT = Mcin , (49) 

and 

3p = MPu . (50) 
The quantities p and T would be the pressure and temperature if the pressure 
tensor and heat flux were, as in classical kinetic theory, connected to the 
fluxes defined with peculiar velocity C, disregarding the convective transport 
involved by the flux J^. 

It can be shown in a similar way to that outlined previously, that p satisfies 

p = nkT . (51) 



15 



In addition, we obtain from equation (126|) 



3p = MF'u - 



1 




(52) 



Mn 



{nvY 



which leads to 



p = p + 




(53) 



and so, from equations (l32l) and (15T1) . 



r 



T + 



P 




(54) 



V 



It appears, therefore, that the pressure and temperature defined within our 
new framework are different from the conventional equilibrium scalars through 
equations fl53|) and fl54l) . The difference may be attributed to the contribution 
of mass-density variations to the pressure and the temperature within the 
fluid, which in turn may be regarded as non-equilibrium contributions to the 
pressure and temperature. The term "non-equilibrium" associated with T' 
and p' refers to the explicit contribution of an additional flux due to density 
variations in their respective expressions, when compared with the expressions 
for T and p. 

An expansion of gas during which the gas is allowed neither to exchange heat 
with its surrounding nor to do external work is called "free expansion", and 
is an irreversible process. As the gas has to expand into its allowed volume, 
work is done internally. Such an internal process within a gas subject to density 
variations is predicted with our new internal energy equation (130!) : while the 
internal energy is constant (i.e. while the non-equilibrium temperature T' 
is a constant), internal work can be done through J„. This work done within 
the gas contributes to changes in the thermal agitation of the particles, i.e. 
changes in T. Moreover, in equations fl53l) and fl54|) the work term due to 
Jy is similar to the classical thermodynamic work [pressure x volume], but 
expressed in terms of density variations, and the work done is transferred to 
the gaseous particle agitation instead of into an external work. So, the new 
contributions of J^, to temperature and pressure resemble a "free expansion" 
heating in which the volume variation is given by the local density gradient. 

If the contribution of J^, is disregarded and the thermal agitation temperature 
of the particles is taken as the thermodynamic parameter and then related 
directly to the internal energy of the particle, as in equation (??) — in other 
words, if we take the classical fluid mechanics viewpoint — then constant 
internal energy implies constant temperature of the gas: the free expansion 
effect should not exist. However, temperature variations due to free expansion 
are well known and experimentally demonstrated. 
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A stationary gas under uniform mechanical (not necessarily free) expansion 
or contraction has Um = 0, so J^, = nvU^. Then equations (l53l) and fl5^ yield 

P' = P+^U^, (55) 

and 

T' = T + ^Ul . (56) 

This means the total pressure of the gas during this expansion or compres- 
sion, when compared with the pressure within the same gas under equilibrium, 
should be augmented by a contribution depending on the square of the ex- 
pansion speed. This additional pressure is a type of dynamic pressure (i.e. the 
pressure given by p[/^/2), where the velocity involved is an expansion velocity 
instead of the mass-velocity. Likewise, the actual temperature of the fluid is 
not simply given by the heat energy corresponding to its internal energy but 
is supplemented by the work done during expansion or contraction. 

It is the case that the differences appearing in equations ( l53l) and ( IMI) between 
the standard and the new definitions of the thermodynamics properties, will 
be negligible for most ordinary flows. 

Two examples will help illustrate the new non-equilibrium thermodynamic 
parameters. Let us first consider a mechanical compression (or an expansion 
against a piston) of argon gas. The molecular mass of argon is M = 6.63 x 
10~^^ kg and the Boltzmann constant is k = 1.38 x 10~^^ J/K. At a piston 
compression speed equivalent to = 50 m/s, the difference in temperature 
due to the work term in equation (!56|) is 4.0 K. 

Let us now consider a shock wave in a rarefied gas flow, within which a rapid 
variation of temperature (at constant pressure) exists over a distance of the 
order of a molecular mean free path. Am, which is of the same order of magni- 
tude as the characteristic length defined by the density gradient. In this case, 
in equation the expansion/compression speed may be approximated by 
II {nv)~^Jy II ~ Xm/Tv where is the characteristic time of diffusion of volume 
(or mass-density). Accordingly, the temperature difference in equation is 




AT = r - T = — — . (57) 



If we compare the diffusion time to the diffusion time for energy, then the 
ratio Xm/'Tv is close to the thermal speed and we have {Xm/TvY ~ 3kT/M. 
Consequently, AT ^ T, i.e. the same order of magnitude as the thermal 
speed temperature (the temperature commonly identified with the flow). This 
difference is similar to some of the temperature differences observed in high 
speed rarefied gas flows during space vehicle re-entry, for which the classical 
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theory of surface temperature jump is known to have serious shortcomings 



5.3 Interpretation of the new convective transport 

The two elementary transport processes admitted so far within fluids are: 
the diffusion attributed to purely random motion of the mass-particles, and 
transport due to the mass-motion or flow. The additional transport process 
due to macroscopic expansion or compression that we posit in this paper could 
be understood as follows. 

To an observer in a fixed reference frame, a fluid could be observed to transport 
energy, mass, or momentum from one point to another by simply undergoing 
expansion or compression. Evidently this cannot be classified as a pure diffu- 
sion process, because the whole fluid should be at rest. It cannot also be clearly 
identified with classical convective transport due to the mass-motion, which is 
physically the flowing of the fluid. This unclassified transport process involves 
the evolution of the domain of the fluid — it is due to the change in volume 
of the fluid. It is a transport process with characteristics of both elementary 
diffusion and convective transport, but distinctively different from each; it is 
similar to what is described in reference [I?] as a "volume transport" . 



6 Conclusions 

In our previous paper we reconsidered the representation of fluids and intro- 
duced a new volume-based kinetic approach for gases. The set of macroscopic 
conservation equations derived from this approach comprises four equations 
rather than the usual three. An evolution equation purely of the mass-density 
is added to the set of three conservation equations. We therefore distinguished 
between conservation of mass, which becomes a conservation of number of 
particles, and the variation of the mass-density as a thermodynamic parame- 
ter to account for any variation in the volume occupied by an element of fluid. 
We also distinguished between compressibility arising from fluid accelerations 
and compressibility arising simply from the physical properties of fluids. 

In this present article we have proposed constitutive relations for the pressure 
tensor and the heat flux to accompany our volume-based macroscopic con- 
servation equations. An unconventional convective form of transport leads us 
to re-classify the diffusive and convective fluxes. This extra convective mode 
of transport can be viewed as a local macroscopic expansion or compression 
motion of the fluid characterized by a flux of volume that classical modelling 
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does not account for. 

As a result, we have obtained a hydro dynamic model in which the contri- 
butions due to a flux of volume appear in the momentum and energy equa- 
tions. Furthermore, the thermodynamic parameters of the fluid (temperature 
and pressure) appear with non-equilibrium contributions depending on density 
variations. 

We propose that our model would show most departure from conventional fluid 
mechanics in micro flow cases, where compressibility effects are significant in 
spite of the low Mach number. It could also be tested on fiows which feature 
large density and temperature variations, such as shock waves. Future work 
should also include: 

• investigating the boundary conditions to accompany our new set of equa- 
tions; 

• analysis of higher-order heat flux and pressure tensor derivations within 
our new kinetic approach — some classical expressions for these (such as 
the Burnett equations) suffer from reference frame dependence, but in our 
model the volume occupied by the molecules is treated separately from the 
fixed reference volume element. This could mitigate some of these problems 
and would, at least, yield interesting new Burnett-like expressions. 
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A The classical view of pressure tensor and heat flux 



In classical kinetic theory, macroscopic expansion or compression is not in- 
volved in the flow description. It is understood that any macroscopic motion 
is defined only by the mass- motion, Um, and so purely random motion would 
be defined by the peculiar velocity C in equation ([8]). In this conceptual frame 
the pressure tensor is given by P and the heat fiux by q [2|3] . Then, in ac- 
cordance with diffusion theory, these quantities are modelled by the following 
standard constitutive relations: 



MPij =p6ij - yU 

M(i=-KhVT 



+ 



dX,, 



' dX, ''^ ' 



(A.l) 



with n the dynamic viscosity, the heat conductivity, rj the bulk viscosity 
[3]. The stress tensor in equations ( lA.ll) is evidently written with the velocity 
Um, as the macroscopic motion associated with the peculiar velocity C is the 
mass velocity. 

In fact, the stress tensor appearing in equations (lA.ip does not account for 
all types of deformations within the fiuid. This is because Um misses any 
expansion/compression of a fiuid element. Strictly speaking, this stress tensor 
accounts only for the "deformation of the mass- velocity" ; it does not account 
for all the deformations of a fiuid element due purely to the changes in the 
fiuid physical properties, e.g. deformation due to heating and due to density 
changes. 

It has long been recognized that the deformation involved in equations ( lA.ll) 
takes place at constant temperature and infers negligible effects from the elas- 
tic nature of fluids [18]. It cannot, therefore, efliciently describe flows in which 
thermo-elastic effects are not negligible. Recently, Brenner used experimental 
results of thermophoretic motion to suggest replacing the velocity in the ex- 
pression of the stress tensor by another velocity that includes these missing 
deformations [T5]. A similar remark can be made for the heat flux expressed 
in equations flA.l|) which also appears to miss effects due to the elastic aspect 
of fluids. 

The parts missing from the current fluid description appear to arise from the 
treatment of elements such as a "fluid particle" and a "control volume of 
fluid" in the derivation of the fluid mechanics set of equations. A fluid particle 
should not be treated as a simple mass-point in which only the mass and the 
mass-velocity are investigated: its volume, and changes in volume, should be 
taken into account in the description. 
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Unlike the classical description, the new stress tensor obtained in our volume- 
based kinetic approach encompasses the elastic nature of the gaseous fluid, 
and allows for thermo-elastic deformations in the fluid. 
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